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NOTE  ON  THE  NUMSilCAL  TREATNEaiT  OF  SECOND-ORDEH  DIFFESIENTUL  EQUATIONS 

by 

RICHARD  UTTQl  AND  HERMAN  KAHN 
(The  RAND  Corporation) 

A  variety  of  aethode  exiat  for  treating  nuaerlcally  aecond- 
order  differential  equatlona  of  the  general  fora 

’  T  («,«*^  ,  (1) 

aupoleaented  by  boundary  conditions 

(l)  0{x)  *  a  at  X  =  a 

=  p  at  X  =  a 
dx 

These  aethoas  have  a  conaon  defect*  however*  being  all  basically 
extranolatlcn  nrocedures  which  Involve  rerlaclng  the  second* order 
differential  by  a  suitable  second  dlfferencec  This  feature  of 
"extranolatloo"  throws  soae  question  on  the  reliability  of  error 
estiaates  deduced  froa  these  aethods*  The  situation  would  be  more 
satisfactory  if  the  solution  of  equation  (l)  could  bo  perforaed  by 
a  coanletely  illicit  scheme*  It  has  been  found  that  such  a  sliarle 
nossiblllty  exists. 

. — ^  The^^reeeat  -scheay  ^roceede  by' ''converting  tiie  differential 
equation  into  an  Integral  equation*  which  Is  effected  by  a  double 


/ 


Interration.r’t*  1 ) r 


r>x 


/( x)  =  a  r  (x-a) 


dy 


dz  F 


(2  »/(*)] 


(2) 


c/ a  Cya 

Th*  tquatlcn  is  net  conveniently  nandltd  In  this  formj  hcwcver» 
an  lrit'‘Lratlon  by  narts  of  the  last  term  fc^lves 

Px 

Ax)  =  a  Ax-a)  \  d>(x-y)  F  [;’,Ay)j  .  (3) 

Ju 

In  this  fox^  thi:  equation  is  particularly  well  suited  for  numerical 
treatrunt.  There  an  two  reasons  for  this:  flr^t,  that  while  Ax) 
at  X  occur-  on  botn  sides  of  the  equation,  it  occurs  with  zero 
wtif  ht  or.  t:.e  ri£,ht  hand  side  and  hence  admits  of  an  explicit  detei^ 
riina'ion  cf  Ax),  and  second,  tne  Integral  in  equation  (3)  Involves 
X  crjy  Iri  as:  Iriessentlal  way  so  that  the  lnteL,ral  niay  be  separated 
into  t-r  intetials,  neither  of  which  conUtiris  x  in  the  intefrand. 
Tnls  means  U.at  equatlcr  (j)  may  h-e  Intet^rated  explicitly  with  a 
slc.plf  accuaulaticn  ever  values  ol  /(x)  on  tne  interval  {a,x). 

In  the  riumericul  treatment  ol  equation  (3)  tne  intcjjral  is 
rf- laced  a  sum.  If  the  sum  is  j  Iver  bj  a  k-pcint  Intet^ratlon 
formula,  tquahlcn  (3)  reduces  to  tiit;  difference  equation 


Ax^)  = 


a 


,  lx  -a) 


X  P(x  )  -  <  (x  ) 
n  n  n  n  1 


(4) 


where 


J=n-k  *“ 


4A 

W(x  )  =  ^  r-,x.  F  x.,/(x.)  ♦  m(x  )  )  (6) 

r  ^ — y  ^  J  J  j  r.-K 

J=tj-k 

and  the  cc€  fflclonts  a.  arise  frcm  the  rpeclftc  Iritei^rsitlon  rcrnula 
which  Is  beint  used.  The  latter  also  datennlnes  the  uccura -v  cf 
tl.e  lr.t«  t  ration.  For  the  tranezoidal  rule,  the  error  is  cf  the  or¬ 
der  cf  the  squart  cf  the  iritex*val  size,  jnlle  for  Slinson's  rule  it 
1e  tm-  feurth  newer,  wiiich  is  equivalent  to  that  cl  the  IferlilaT 
Runi;t— Kutta  method. 

In  erdfc-r  to  start  the  Intetraticn  of  the  al'lerence  equation  (..) 
it  le  necessarj'  U-  a«-iecif>  the  vaiuni  of  /(x),  rCx)  and  ^(x)  at  ti." 
first  (k-i.)  Dclr.tB  cf  the  lr.tot;ratlon.  Tnls  requir*  a  an  anuljtl  : 
armroxiaatlon  about  the  st^rtintj  nclrit  x  =  a.  This  Is  accnnilshed 
by  first  eXTiaridlnr  F  [x,/(x)]  as  a  newer  series  in  x  and  Uien  dttei- 
iulnini,  /(x),  P(x),  and  WCx)  frem  tnls  serifs.  Th.us  let 


F  (x,/(x)]  =  /_,  77 


wilt  re 


A.  =  —  f  lx,/{x)] 

QX^ 


Ttie  evaluation  oC  ia  carried  cut  dlre<  with  the.  aid  cl  U  e 


expn^aaicns 


/(a)  =  a 
/•(a)  =  r 
/"(a)  =  F(a,a) 
/••(x)  =  F  (x,/(x)] 


(9) 

(10) 

(11) 

(12) 


which  rdlrw  limnedlatcly  from  equation  (l)  and  Ita  associated 
Ix^ur.tiari  ccndl lions#  equation  (?'  when  substituted  Irtc 


F(x)  =  \  dyF  [y,/(j)] 

-'a 

X 

qCx)  =  V  dyj-  F  [y, /(>•)] 


(13) 


(U) 


and 

Ax)  =  a  ♦  , :(x-a)  fx?(x)  -  t«(x)]  ,  (1^.) 

ccnctitute  ti.e  required  frnniilas  to  Initiate  the  intet'ratlct.,  k'f 
coiL'sc,  one  retains  only  these  letus  cf  eq’watlcn  (7)  which  are  con¬ 
sistent  will,  trie  accuracy  cf  trie  Iritefratlcn  fcrmula. 

The  metrod  presented  here  has  been  a^nlied  In  detalJ  tc  the 
sneciflc  case  In  widch  (l)  Is  ths  well  knewn  Fermi-Tliomas  eqoatlcn. 
In  addl’ion,  fer  Inis  casey  the  ar.alytic  sim»^llclty  of  trie  intecraJ 
crmulatlon  •as  been  exolcited  to  obtain  anoroxiaate  solutions. 


Tnese  results  will  be  re->orted  elsewhere. 


